Finite Size Scaling for Criticality of the Schrodinger Equation 
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By solving the Schrodinger equation one obtains the whole energy spectrum, both the bound and 
the continuum states. If the Hamiltonian depends on a set of parameters, these could be tuned to 
a transition from bound to continuum states. The behavior of systems near the threshold, which 
separates bound-states from continuum states, is important in the study of such phenomenon as: 
ionization of atoms and molecules, molecule dissociation, scattering collisions and stability of matter. 
In general, the energy is non-analytic as a function of the Hamiltonian parameters or a bound-state 
does not exist at the threshold energy. The overall goal of this chapter is to show how one can 
predict, generate and identify new class of stable quantum systems using large-dimensional models 
and the finite size scaling approach. Within this approach, the finite size corresponds not to the 
spatial dimension but to the number of elements in a complete basis set used to expand the exact 
eigenfunction of a given Hamiltonian. This method is efficient and very accurate for estimating the 
critical parameters, {Xi}, for stability of a given Hamiltonian, H(\i). We present two methods of 
obtaining critical parameters using finite size scaling for a given quantum Hamiltonian: The finite 
element method and the basis set expansion method. The long term goal of developing finite size 
£1^ scaling is treating criticality from first principles at quantum phase transitions. In the last decade 

considerable attention has concentrated on a new class of phase transitions, transitions which occur 
at the absolute zero of temperature. These are quantum phase transitions which are driven by 
quantum fluctuations as a consequence of Heisenberg's uncertainty principle. These new transitions 
are tuned by parameters in the Hamiltonian. Finite size scaling might be useful in predicting the 
quantum critical parameters for systems going through quantum phase transitions. 



> 

m 

ON 

o 
o 



x 



kais@purdue.edu 



I. INTRODUCTION 



Weakly bound states represent an interesting field of research in atomic and molecular physics. The behavior 
of systems near the binding threshold is important in the study of ionization of atoms and molecules, molecule 
dissociation and scattering collisions. Moreover, the stability of atomic and molecular systems in external electric, 
magnetic and laser fields is of fundamental importance in atomic and molecular physics and has attracted considerable 
experimental and theoretical attention over the past decades [US]. A superintense laser field can change the nature 
of atomic and molecular systems and their anions; the stabilization in superstrong fields is accompanied by splitting 
of the electron distribution into distinct lobes, with locations governed by the quiver amplitude and polarization of 
the laser field. This localization markedly alters electron-nucleus interactions as well as reduces electron-electron 
repulsions and hence suppresses autoionization. In molecules, it can also enhance chemical bonding. This localization 
markedly reduces the ionization probability and can enhance chemical bonding when the laser strength becomes 
sufficiently strong and can give rise to new stable multiply charged negative ions such as H , He - and HJT [31 [SHI] . 

In general, the energy is non-analytical, an analytic function is a function that is locally given by a convergent 
power series, as a function of the Hamiltonian parameters or a bound-state does not exist at the threshold energy. 
It has been suggested for some time, based on large-dimensional models, that there are possible analogies between 
critical phenomena and singularities of the energy [8TfT0] . 

Phase transitions are associated with singularities of the free energy. These singularities occur only in the thermo- 
dynamic limit [111 I12j where the dimension of the system approaches infinity. However calculations are done only on 
finite systems. A Finite Size Scaling (FSS) approach is needed in order to extrapolate results from finite systems to 
the thermodynamic limit |13|. FSS is not only a formal way to understand the asymptotic behavior of a system when 
the size tends to infinity, but a theory that also gives us numerical methods |TH - f2"D] capable of obtaining accurate 
results for infinite systems by studying the corresponding small systems 21 32 . Applications include expansion in 
Slater-type basis functions [30 , Gaussian-type basis functions [33 j and recently finite elements |34j. 

II. CRITIC ALITY FOR LARGE-DIMENSIONAL MODELS 

Large dimension models were originally developed for specific theories in the fields of nuclear physics, critical 
phenomena and particle physics[351 135]. Subsequently, with the pioneering work of HerschbachjTni |37], they found 
wide use in the field of atomic and molecular physics 38 1. In this method one takes the dimension of space, D, as 
a variable, solves the problem at some dimension fl / 3 where the physics becomes much simpler, and then uses 
perturbation theory or other techniques to obtain an approximate result for D = 3|10). 

It is possible to describe stability and symmetry breaking of electronic structure configurations of atoms and 
molecules as phase transitions and critical phenomena. This analogy was revealed by using dimensional scaling 
method and the large dimensional limit model of electronic structure configurations [59TH2"] . 

To study the behavior of a given system near the critical point, one has to rely on model calculations which are 
simple, capture the main physics of the problem and which belong to the same universality class |161 117j . Here we 
will illustrate the phase transitions and symmetry breaking using the large dimension model. In the application of 
dimensional scaling to electronic structure, the large-D limit reduces to a semi-classical electrostatic problem in which 
the electrons are assumed to have fixed positions relative both to the nuclei and to each other in the D-scaled space 10J . 
This configuration corresponds to the minimum of an effective potential which includes Coulomb interactions as well 
as centrifugal terms arising from the generalized D-dependence kinetic energy. Typically, in the large-D regime the 
electronic structure configuration undergoes symmetry breaking for certain ranges of nuclear charges or molecular 
geometries [3B], 

In order to illustrate the analogy between symmetry breaking and phase transitions we present as an example: 
the results for the two-electron atoms in the Hartree-Fock (HF) approximation [39]. In the HF approximation at the 
D — > oo limit, the dimensional-scaled effective Hamiltonian for the two-electron atom in an external weak electric 
field £ can be written as [171H5] . 
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where r% and are the electron-nucleus radii, and Z is the nuclear charge. The ground state energy at the large-D 
limit is then given by E ao (Z,£) = min/ ri)r2 \ Hoc- 

In the absence of an external electric field, £ — 0, Herschbach and coworkers [15] have found that these equations 
have a symmetric solution with the two electrons equidistant from the nucleus, with r% = r-i = r. This symmetric so- 
lution represents a minimum in the region where all the eigenvalues of the Hessian matrix are positive, Z > Z c = 
For values of Z smaller than Z c , the solutions become unsymmetrical with one electron much closer to the nucleus 
than the other {r\ 7^ r%). In order to describe this symmetry breaking, it is convenient to introduce new variables 
(r, r\) of the form: T\ = r; r-i = (1 — rf)r, where r) = (n — T2)fr\ 7^ measures the deviation from the symmetric 
solution. 

By studying the eigenvalues of the Hessian matrix, one finds that the solution is a minimum of the effective potential 
for the range, 1 < Z < Z c . We now turn to the question of how to describe the system near the critical point. To 
answer this question, a complete mapping between this problem and critical phenomena in statistical mechanics is 
readily feasible with the following analogies: 

• nuclear charge (Z) o temperature (T) 

• external electric field (£) -H* ordering field (h) 

• ground state energy (E OQ (Z,£)) «-» free energy (f(T,h)) 

• asymmetry parameter (77) O order parameter (m) 

• stability limit point (Z c ,£ = 0) ■<-» critical point (T c , h = 0) 

Using the above scheme, we can define the critical exponents (f3,a,6 and 7) for the electronic structure of the two 
electron atom in the following way: 



r,(Z,£ = 0) 
£oo(Z,£ = 0) 

as l£=o 
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I AZ |" 
V S sgn(ri) 
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AZ 0- 
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T) -> 

AZ — > 



(2) 



where AZ = Z — Z c . These critical exponents describe the nature of the singularities in the above quantities at 
the critical charge Z c . The values obtained for these critical exponents are known as classical or mean- field critical 
exponents: /3 = \ ; a = 2, ; 5 = 3; 7 = 1. 



This analogy between symmetry breaking and phase transitions was also generalized to include the large dimen- 
sional model of the N-electron atoms [40], simple diatomic molecules [41 [ 143] . both linear and planar one-electron 
systems [H] as well as three-body Coulomb systems of the general form Ai?^4|^4"]. 

The above simple large-D picture helps to establish a connection to phase transitions. However, the next question 
to be addressed is: How to carry out such an analogy to D — 3?. This question will be examined in the subsequent 
sections using the finite size scaling approach. 



III. FINITE SIZE SCALING: A BRIEF HISTORY 



Ice tea, boiling water and other aspects of two-phase coexistence are familiar features of daily life. Yet phase 
transitions do not exist at all in finite systems! They appear in the thermodynamic limit: The volume V — > 00 and 
particle number N — > 00 in such a way that their ratio, which is the density p = N/V , approaches a finite quantity. In 
statistical mechanics, the existence of phase transitions is associated with singularities of the free energy per particle 
in some region of the thermodynamic space. These singularities occur only in the thermodynamic limit^Tll Q2] . This 
fact could be understood by examining the partition function Z. 

Z= J2 e- E ^/ k * T , (3) 

microstate O 

where E(£l) is the energies of the states, fee is the Boltzmann constant and T is the temperature. For a finite 
system, the partition function is a finite sum of analytical terms, and therefore it is itself an analytical function. The 
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Boltzmann factor is an analytical function of T except at T = 0. For T > 0, it is necessary to take an infinite number 
of terms in order to obtain a singularity in the thermodynamic limit|ll | 112 ) . 

In practice, real systems have a large but finite volume and particle numbers (N ~ 10 23 ), and phase transitions 
are observed. More dramatic even is the case of numerical simulations, where sometimes systems with only a few 
number (hundreds, or even tens) of particles are studied, and "critical" phenomena are still present. Finite size 
scaling theory, which was pioneered by Fisher [IB], addresses the question of why finite systems apparently describe 
phase transitions and what is the relation of this phenomena with the true phase transitions in corresponding infinite 
systems. Moreover, finite-size scaling is not only a formal way to understand the asymptotic behavior of a system 
when the size tends to infinity. In fact, the theory gives us numerical methods capable of obtaining accurate results 
for infinite systems by studying the corresponding small systems (see [15H17) and references therein). 

In order to understand the main idea of finite size scaling, let us consider a system defined in a _D-dimensional vol- 
ume V of a linear dimension L (V = L D ). In a finite size system, If quantum effects are not taken into consideration, 
there are in principle three length scales: The finite geometry characteristic size L, the correlation length £, which 
may be defined as the length scale covering the exponential decay e~ r /£ with distance r of the correlation function, 
and the microscopic length a which governs the range of the interaction. Thermodynamic quantities thus may depend 
on the dimensionless ratios £/a and L/a. The finite size scaling hypothesis assumes that, close to the critical point, 
the microscopic length drops out. 

If in the thermodynamic limit, L — > oo, we consider that there is only one parameter (say temperature T) in the 
problem and the infinite system has a second order phase transition at a critical temperature T c , a thermodynamic 
quantity G develops a singularity as a function of the temperature T in the form: 



G(T) = lim G L (T) ~ \T - T c \~ p , (4) 

whereas it is regular in the finite system, Gl(T) has no singularity. 

When the size L increases, the singularity of G(T) starts to develop. For example, if the correlation length diverges 

clt Tq cLS I 



£(T)= lim UT)^\T-T C \~ U , (5) 

L— >oo 

then £,l{T) has a maximum which becomes sharper and sharper, then FSS ansatz assumes the existence of scaling 
function Fk such that: 



G L {T) ~ G{T)F K (^^) , (6) 

where Fx{y) ~ y p ^ v for y ~ + . Since the FSS ansatz, Eq. should be valid for any quantity which exhibits an 
algebraic singularity in the bulk, we can apply it to the correlation length £ itself. Thus the correlation length in a 
finite system should have the form: 



i L {T)^L^{L^\T-T c \). (7) 

The special significance of this result was first realized by Nightingale [51] , who showed how it could be reinterpreted as 
a renormalization group transformation of the infinite system. The phenomenological renormalization (PR) equation 
for finite systems of sizes L and L' is given by: 



and has a fixed point at T^ L ' L \ It is expected that the succession of points |t( L ' L H will converge to the true T c ir 



the infinite size limit. 
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The finite-size scaling theory combined with transfer matrix calculations had been, since the development of the 
phenomenological renormalization in 1976 by Nightingale [51], one of the most powerful tools to study critical phe- 
nomena in two-dimensional lattice models. For these models the partition function and all the physical quantities of 
the system (free energy, correlation length, response functions, etc) can be written as a function of the eigenvalues of 
the transfer matrix [5^]. In particular, the free energy takes the form: 

f(T) = -Tin Ai (9) 

and the correlation length is: 

£C0 = , , * s (io) 

ln(A 2 /Ai) 

where Ai and A2 are the largest and the second largest eigenvalues of the transfer matrix. In this context, critical 
points are related with the degeneracy of these eigenvalues. For finite transfer matrix the largest eigenvalue is isolated 
(non degenerated) and phase transitions can occur only in the limit L — ¥ 00 where the size of the transfer matrix 
goes to infinity and the largest eigenvalues can be degenerated. In the next section, we will see that these ideas of 
finite size scaling can be generalize to quantum mechanics, in particular addressing the criticality of the Schrodinger 
equation. 

IV. FINITE SIZE SCALING FOR THE SCHRODINGER EQUATION 

The finite size scaling method is a systematic way to extract the critical behavior of an infinite system from analysis 
on finite systems [30]. It is efficient and accurate for the calculation of critical parameters of the Schrodinger equation. 
Let's assume we have the following Hamiltonian: 

H = Ho + V x (11) 

where Hq is A-independent and V\ is the A-dependent term. We are interested in the study of how the different 
properties of the system change when the value of A varies. A critical point, A c , will be defined as a point for which 
a bound state becomes absorbed or degenerate with a continuum. 



Without loss of generality, we will assume that the Hamiltonian, Eq. (11 1, has a bound state, E\, for A > A c which 
becomes equal to zero at A = A c . As in statistical mechanics, we can define some critical exponents related to the 
asymptotic behavior of different quantities near the critical point. In particular, for the energy we can define the 
critical exponent a as: 

E x „ (A-A C ) Q . (12) 

A ->• A+ 

The existence or absence of a bound state at the critical point is related to the type of the singularity in the 
energy. Using statistical mechanics terminology, we can associate "first order phase transitions" with the existence 
of a normalizable eigenfunction at the critical point. The absence of such a function could be related to "continuous 
phase transitions" [3TJJ . 

In quantum calculations, the variational method is widely used to approximate the solution of the Schrodinger 
equation. To obtain exact results one should expand the exact wave function in a complete basis set and take the 
number of basis functions to infinity. In practice, one truncates this expansion at some order N. In the present 
approach, the finite size corresponds not to the spatial dimension, as in statistical mechanics, but to the number of 
elements in a complete basis set used to expand the exact eigenfunction of a given Hamiltonian. We will compare 
two methods to obtain the matrix elements needed to apply the FSS ansatz. The size of our system for the basis set 
expansion will correspond to the dimension of the Hilbert space. For a given complete basis set <&„, the ground-state 
eigenfunction has the following expansion: 



*A=^fln(A)^„, (13) 



G 



where n is the set of quantum numbers. We have to truncate the series at order N and the expectation value of any 
general operator O at order N is given by: 



N 



(Of = 5>WaWO n>ro , (14) 



where O n . m are the matrix elements of O in the basis set {ip n }- 

For the finite element method (FEM), the wavefunction ip n (r) in the n-th element is expressed in terms of local 
shape functions. For our calculations, we use Hermite interpolation polynomials with two nodes and three degrees of 
freedom. This choice ensures the continuity of the wavefunction and its first two derivatives. Then in n-th element 
the wavefunction is 134 : 



Mr) = £ [MrW n + MrW': + U^'n] , (15) 



with a indicating the nodal index of the element; i = 1 for the left and i — 2 for the right border of the element. The 
functions </>i(r), 4>i(r), and 4>i{r) are the (fifth degree) Hermite interpolation polynomials. Then ip^, ip^, and ip n l are 
the undetermined values values of the wavefunction and its first and second derivative on the nodal points. The size 
for the case of solving the equation with the FEM will be the number of elements used. 

Since (0) A is not analytical at A = A c , then we define a critical exponent, fiQ, if the general operator has the 
following relation: 

(0) x ^(X-X c r° for A^A+, (16) 

where A — 5- A^j" represents taking the limit of A approaching the critical point from larger values of A. As in the 
FSS ansatz in statistical mechanics [16| [53| , we will assume that there exists a scaling function for the truncated 
magnitudes such that: 

(0)W~{0) x F (N\\-X c n (17) 
with the scaling function Fq being particular for different operators but all having the same unique scaling exponent 

To obtain the critical parameters, we define the following function: 

ln((Q)f/(Q)f) 

A o(X,N,N)- Hn//n) . (18) 



At the critical point, the expectation value is related to N as a power law, (O) ~ N^ ^, and Eq. (18) becomes 
independent of N. For the energy operator O = H and using the critical exponent a for the corresponding exponent 
fio w e have: 

A H (\ C ;N,N') = ^-. (19) 
In order to obtain the critical exponent a from numerical calculations, it is convenient to define a new function|30j: 

MW) = — f „ (20) 

which at the critical point is independent of N and N' and takes the value of a. Namely, for A = A c and any values 
of N and TV' we have 

T a (\ Cl N 1 N') = a, (21) 
and the critical exponent v is readily given by Eq. (|19|. 
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V. THE HULTHEN POTENTIAL 



To illustrate the application of the FSS method in quantum mechanics, let us give an example of the criticality 
of the Hulthen potential. The Hulthen potential behaves like a Coulomb potential for small distances whereas for 
large distances it decreases exponentially so that the "capacity" for bound states is smaller than that of Coulomb 
potential. Thus, they have the same singularity but shifted energy levels. They always lie lower in the Coulomb case 
than in the Hulthen case, where there remains only space for a finite number of bound states [55]. Here, we present the 
FSS calculations using two methods: finite elements and basis set expansion; each used to obtaining quantum critical 
parameters for the Hulthen Hamiltonian. First, we give the analytical solution, then FSS with basis set expansion 
and finite element solution. 



A. Analytical Solution 

The Hulthen potential has the following form |541 155] : 

V U = -*T=T* (22) 

where A is the coupling constant, and a is the scaling parameter. For small values of r/a the potential V(r) — > — -X/r, 
whereas for large values of r/a the potential approaches zero exponentially fast, therefore the scale a in the potential 
regulates the infinite number of levels that would otherwise appear with a large-distance Coulomb behavior. 

Shrodinger radial differential equation in the dimensionless variable r — r/a becomes: 

Hp + c- ,+ *i^>*-°- < 23 > 

We only consider the case for I — for the Hulthen potential. Here we used the abbreviations a 2 = —Ea 2 > 
(in atomic units m = Ti — 1). The complete solutions for the wavefunctions are written in term of hypergeometric 
functions as follows (55 : 



x = N e- ar (l - e- r ) 2 F 1 (2a + 1 + n, 1 - n, 2a + 1; e~ r ), (24) 

where the normalization factor is given by No = [a(a + n)(2a + n)] i [T(2a + n)/T(2a + l)T(n)]. It follows that the 
energy levels are given by: 

_ 1 (2A-n 2 ) 2 , 

We can make the following comments concerning the energy levels obtained for the Hulthen potential. There exists 
a critical value for the coupling A c to have the given energy levels, A c = n 2 /2. It follows directly from the first 
observation that the number of levels n max allowed is finite and it depends on the size of the coupling constant 
n max — 2A. As A — > oo the potential is well behaved, which can be seen as follows: In this limit we get the obvious 
inequality a 2 <C 2A =>• v2A ~ n. It follow that we can set a w in Eq. (24 1 to obtain: 

Xc^o = (1 - e- r ) 2 F!(l + n, 1 - n, 1; e~ r ), (26) 

which is the wave function at threshold. This wave function is not normalizable as expected when the energy exponent 
a = 2, E ~ (A — \ c ) a . For the ground state, the asymptotic limit of the probability density for r >> 1 and A — s- A c 
becomes: 

P(r)~e"^, £~|A-A c r, (27) 



with a characteristic length £ and exponent v = 1. The Hulthen potential has a finite capacity determined by the 
critical coupling, A c . The potential admits bound states between the range of values for the coupling: A = [1/2, oo). 
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FIG. 1. Plot of r a , obtained by FSS method, as a function of A. Using the number of basis N from 8 to 48 in steps of 2. For 
FEM the number elements used were from 100 to 380 in steps of 20. 



B. Basis Set Expansion 

For the Hulthen potential, the wavefunction can be expanded in the following Slater basis ( see Chapter 7 for details 

*»(r) = J , Ji{t+» e ~ r/ ' L »W- (28) 
y (n + l)(n + 2) 

(21 

Ln (r) is the Laguerre polynomial of degree n and order 2. The kinetic term can be obtained analytically. However, 
the potential term need to be calculated numerically |5 7j . 

Figure 1, show the results for the plot r Q (A, N, N') as a function of A with different N and N' , all the curves will 
cross exactly at the critical point. 



C. Finite Element Method 

The FEM is a numerical technique which gives approximate solutions to differential equations. In the case of 
quantum mechanics, the differential equation is formulated as a boundary value problem [6TJ [62]. For our purposes, 
we are interested in solving the time-independent Shrodinger equation with finite elements. We will require our 
boundary conditions to be restricted to the Dirichlet type. For this problem, we will use two interpolation methods: 
linear interpolation and Hermite Interpolation polynomials to solve for this potential. 

We start by integration by parts and impose the boundary conditions for the kinetic energy and reduce it to the 
weak form|34j: 

- I r 2 ip*' (r)ip' (r)dr. 
2 Jo 

(29) 

For the potential energy: 

poo —e_r r 

I r 2 ip*(r)ip(r)X- —dr. 

Jo 1 - e 

(30) 

We calculated the local matrix elements of the potential energy by using a four point Gaussian Quadrature to 
evaluate the integral. We set the cutoff for the integration to r c . To include the integration to infinity, we added an 
infinite element approximation. To do so, we approximate the solution of the wave function in the region of [r c , oo) 
to be an exponentially decaying function with the form ip( r ) = ^( r c) &~ r ■ 

The local matrices are then assembled to form the complete solution and by invoking the variational principle on 




9 



the nodal values ipi we obtain a generalized eigenvalue problem representing the initial Schrodinger equation: 

H lJ \^)=eU lJ \^ J ). (31) 



Solution of Eq. (31 ) is achieved using standard numerical methods (see Chapter 10 for details 
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FIG. 2. Extrapolated values for the critical exponents and the critical parameter A. The solid read dots at 1/N — are 
the extrapolated critical values. The left side is the basis set method while the right is the FEM with Hermite interpolation 
polynomials. 



D. Finite Size Scaling Results 




FIG. 3. Data collapse study of the basis set method and FEM. The left is the basis set method and the right being the FEM. 



The finite size scaling equations are valid only as asymptotic expressions, but unique values of A c , a, and v can be 
obtained as a succession of values as a function of N. The lengths of the elements are set h = 0.5. The plots of T a , 
figure [T] the basis set expansion is giving values very close to the analytical solution of the Hulthen potential. For the 
plot of T a for the FEM estimation of A c is producing results very close to the exact values using Hermite interpolation. 
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The intersection of these curves indicate the A c on the abscissa. The ordinate gives the critical exponents a (in T a 
plots). In Figure]^ we observed the behavior of the pseudocritical parameters, \i N \ ai N \ vi N \ as a function of 1/N. 
The three curves monotonically converge to limiting values for the Hermite interpolation and the basis set expansion. 

To check the validity of our finite size scaling assumptions, we performed a data collapse[51] calculation of the 
Hulthen potential. In the data collapse analysis, we examine the main assumption we have made in Eq. (17) for the 
existence of a scaling function for each truncated mag nitude (0)^° with a unique scaling exponent v. 

Since the (O)^ is analytical in A, then from Eq. (17) the asymptotic behavior of the scaling function must have 
the form: 

F (x) ~ x-»°/ v . (32) 

For large values of N, at the A c , we have 

(0) {N) (Ac)~ N-^l". (33) 
Because the same argument of regularity holds for the derivatives of the truncated expectation values, we have: 



<9A" 



j\j—(fJ,o-m)/u 



(34) 



A=A C 



(Oy N ' is analytical in A, then using Eq. (34), the Taylor expansion could be written as 



(0) {N) (A) ~ N-^^GoiN^^X - Ac)), (35) 
where Go is an analytical function of its argument. This equivalent expression for the scaling of a given expectation 



value has a correct form to study the data collapse in order to test FSS hypothesis. If the scaling Eq. (17) or Eq. (35) 
holds, then near the critical point the physical quantities will collapse to a single universal curve when plotted in the 
appropriate form (Oy N ^ N^°l v against N 1 ^(X — A c ). If the operator O is the Hamiltonian then we will have a data 
collapse when plotting EaN~ a / v against A^ 1 / l/ (A — A c ). In Figure [SJwe plot the results corresponding to the basis set 
method (right panel) and Hermite interpolation (left panel), which have been calculated with A c = 0.49999, a = 1.9960 
and v = 0.99910 for the basis set method and for the Hermite interpolation we have A c = 0.50000, a = 2.00011 and 
v = 1.000322. The data collapse study do in fact support our FSS assumptions. We have conveniently summarized 
our results for the critical parameters for the analytical, linear interpolation, Hermite interpolation and the basis set 
expansion in table |T] 



TABLE I. Critical Parameters for the Hulthen Potential 
Analytical Linear Hermite Basis Set 
A 0.5 (exact) 0.50184 0.50000 0.49999 
a 2 (exact) 1.99993 2.00011 1.9960 
v 1 (exact) 1.00079 1.00032 0.99910 



We have successfully obtained the critical exponents and the critical parameter for the Hulthen potential using FSS 
with the basis set method and the FEM. The results are in excellent agreement with the analytical solution even for 
the very simplistic linear interpolation used for the FEM calculations. However, the ability of the FEM to describe 
the wavefunction locally in terms of elements affords a very natural way to extend its use for FSS purposes. 
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VI. FINITE SIZE SCALING AND CRITIC ALITY OF M-ELECTRON ATOMS 



Let us examine the criticality of the N-electrons atomic Hamiltonian as a function of the nuclear charge Z . The 
scaled Hamiltonian takes the form: 



M 



1. 



M 



(36) 



i<3 = l 



where are the interelectron distances, and A = 1/Z is the inverse of the nuclear charge. For this Hamiltonian, a 
critical point means the value of the parameter, A c , for which a bound state energy becomes absorbed or degenerate 
with the continuum. 

To carry out the FSS procedure, one has to choose a convenient basis set to obtain the two lowest eigenvalues and 
eigenvectors of the finite Hamiltonian matrix. For M — 2, one can choose the following basis set functions: 



*«M(*i»2a) = ^= (r\r{e-^ +s ^+ r^e"^"^) r k 12 F/(0 ia ,n) (37) 

where 7 and 5 are fixed parameters, we have found numerically that 7 = 2 and 8 = 0.15 is a good choice for the ground 
state[5T], ri2 is the interelectronic distance and Fg(0 12, ft) is a suitable function of the angle between the positions 
of the two electrons #12 and the Euler angles ft = (O,^, 4 r ). This function i*i is different for each orbital-block of 
the Hamiltonian. For the ground state Fo(6i2,CI) = 1 and Fi(9i2,tt) = sin(#i2) cos(O) for the 2p 2 3 P state. These 
basis sets are complete for each ^-subspace. The complete wave function is then a linear combination of these terms 
multiplied by variational coefficients determined by matrix diagonalization [3T]. In the truncated basis set at order N, 
all terms are included such that N > i+ j + k. Using FSS calculations with N = 6, 7, 8, . . . , 13 gives the extrapolated 
values of A c = 1.0976 ± 0.0004 which is in excellent agreement with the best estimate of A c = 1.09766079 using 
large-order perturbation calculations [65 . Since the critical charge Z c = 1/A C ~ 0.91 indicates that the hydrogen 
anion H~ is stable, Z = 1 > Z c . 

For three-electron atoms, M = 3, one can repeat the FSS procedure with the following Hyllerass-type basis set|22j: 

9 ij u m n@i,M)=CA(ri4r$r l ia r&r3 1 e -«Cn+r*) e -*a Xl ) , (38) 

where the variational parameters, a = 0.9 and f3 — 0.1, were chosen to obtain accurate results near the critical charge 
Z ~ 2, xi is the s P m function with spin angular moment 1/2: 



Xi = a(l)0(2)a(3) - /3(l)a(2)a(3), (39) 

C is a normalization constant and A is the usual three-particle antisymmetrizer operator [22] . The FSS calculations 
gives A c = 0.48 ± 0.03. Since Z c ~ 2.08 the anions He - and H are unstable. 

One can extend this analysis and calculate the critical charges for M-electron atoms in order to perform a systematic 
check of the stability of atomic dianions. In order to have a stable doubly negatively charged atomic ion one should 
require the surcharge, S e (N) = N — Z C (N) > 2. We have found that the surcharge never exceeds two. The maximal 
surcharge, S e (86) — 1.48, is found for the closed-shell configuration of element Rn and can be related to the peak of 
electron affinity of the clement N = 85. The FSS numerical results for M-electron atoms show that at most, only 
one electron can be added to a free atom in the gas phase. The second extra electron is not bound by singly charged 
negative ion because the combined action of the repulsive potential surrounding the isolated negative ion and the 
Pauli exclusion principle. However, doubly charged atomic negative ions might exist in a strong magnetic field of the 
order few atomic units, where la.u. = 2.3505 10 9 G and superintese laser fields. 



VII. CONCLUSIONS 



In this chapter, we show how the finite size scaling ansatz can be combined with the variational method to extract 
information about critical behavior of quantum Hamiltonians. This approach is based on taking the number of 
elements in a complete basis set or the finite element method as the size of the system. As in statistical mechanics, 
finite size scaling can then be used directly to the Schrodinger equation. This approach is general and gives very 
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accurate results for the critical parameters, for which the bound state energy becomes absorbed or degenerate with 
a continuum. To illustrate the applications in quantum calculations, we present detailed calculations for the simple 
case of Hulthcn potential and few electron atoms. For atomic systems we have shown that finite size scaling can be 
used to explain and predict the stability of atomic anions: At most, only one electron can be added to a free atom in 
the gas phase. 

Recently, there has been an ongoing experimental and theoretical search for doubly charged negative molecular 
dianionspQ. In contrast to atoms, large molecular systems can hold many extra electrons because the extra electrons 
can stay well separated. However, such systems are challenging from both theoretical and experimental points of 
view. The present finite size scaling approach might be useful in predicting the general stability of molecular dianions. 

The approach can be generalize to complex systems by calculating the matrix elements needed for FSS analy- 
sis by ab initio, density functional methods, orbital free density functional (OF-DFT) [SSI HZ] approach, density 
matrices [68] [69] and other electronic structure methods[70]. The implementation should be straightforward. We need 
to obtain the matrix elements to calculate r a as a function of the number of elements used in solving for the system. 
In the finite element using mean field equations (like Hartree-Fock or Kohn Sham methods) the solution region will 
be discretized into elements composed of tetrahedrons. 

The field of quantum critical phenomena in atomic and molecular physics is still in its infancy and there are many 
open questions about the interpretations of the results including whether or not these quantum phase transitions 
really do exist. The possibility of exploring these phenomena experimentally in the field of quantum dots [71 j and 
systems in superintense laser fields [72"] offers an exciting challenge for future research. This finite size scaling ap- 
proach is general and might provide a powerful way in determining critical parameters for the stability of atomic and 
molecular systems in external fields, and for design and control electronic properties of materials using artificial atoms. 

The critical exponents calculated with finite size scaling indicate the nature of the transitions from bound to 
continuum states. Study of the analytical behavior of the energy near the critical point show that the open shell 
system, such as the lithium like atoms, is completely different from that of a closed shell system, such as the helium 
like atoms. The transition in the closed shell systems from a bound state to a continuum resemble a "first-order phase 
transition", E ~ (A — Ac) 1 , while for the open shell system, the transition of the valence electron to the continuum 
is a "continuous phase transition", E ~ (A — A c ) 2 . For closed shell systems, one can show that H(A C ) has a square- 
integrable eigenfunction corresponding to a threshold energy, the existence of a bound state at the critical coupling 
constant A c implies that for A < A c , E(X) approaches E(X C ) linearly in (A — A c ) as A — > X~ . However, for open shell 
systems, the wave function is not square-integrable at at A c . This difference in critical exponents might be helpful in 
developing a new atomic classification schemes based on the type of phase transitions and criticality of the system. 
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